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Abstract 

The lattice dynamics in Z'^ , d > 1 , is considered. The initial data are supposed 
to be random function. We introduce the family of initial measures > 0} de- 

pending on a small scaling parameter e . We assume that the measures /ig are locally 
homogeneous for space translations of order much less than and nonhomogeneous 
for translations of order . Moreover, the covariance of Hq decreases with distance 
uniformly in e . Given rGM\0, rGR'^,and k>0, we consider the distribu- 
tions of random solution in the time moments t = t je'^ and at lattice points close to 
[r/e] G . The main goil is to study the asymptotics of these distributions as e — )• 
and derive the limit hydrodynamic equations of the Euler or Navier-Stokes type. The 
similar results are obtained for lattice dynamics in the half-space Z^J. . 
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1 Introduction 



One of the central problems in nonequilibrium statistical physics is the derivation of hydrody- 
namic equations of fluid from the microscopic Hamilton dynamics. The main fluid equations 
are the Euler and Navier-Stokes equations. The idea that the Euler equation of fluid dy- 
namics could be derived from microscopic dynamics goes back to Morrey [22]. A systematic 
explanation of some basic ideas and first results are presented in the survey papers by De 
Masi, laniro, Pellegrinotti, and Presutti [3J, by Dobrushin, Sinai and Sukhov [5j, and by 
Spohn [23J. 

One approach is to derive the Euler and Navier-Stokes equations from the Boltzmann 
equation (see, for example, However, the Boltzmann equation is not a microscopic 

model, it should itself be derived as a scaling limit of a more basic model. An alternate 
approach is to study the hydrodynamic behaviour of some simplified or idealized models of 
interacting particles. For models of stochastic dynamics, the results were obtained by Yau et 
al. (see, for example, [IHl [201 EI] , the survey paper by Fritz [12] and the bibliography there). 
Deterministic models, where only initial realizations can be random, are more difficult to 
study in a hydrodynamic framework. The first results were obtained by Dobrushin et al. 
for the one-dimensional hard rods [Ij, and for the one-dimensional oscillators on the lattice 
The main purpose of these models is to show how hydrodynamic behavior arises. 

The present work continues the papers [TH [T7] , where as the model the harmonic crystals 
in are considered. In the harmonic approximation, the crystal is characterized by the 
displacements u{z,t) G M" , z & Z'^ , of the crystal atoms from their equilibrium positions. 
The field u{z,t) is governed by a discrete wave equation. The harmonic crystals can be 
considered as an extension of the model of the infinite chain of one-dimensional harmonic 
oscillators to the many dimensional case. 

The derivation of hydrodynamic equations is connected with the problem of convergence 
to an equilibrium measure. For harmonic crystals, such convergence was proved in [Tl] IT2 | 115] . 
We outline this result. The initial data are assumed to be random function with a distribution 
/io • We suppose that the measure fiQ has zero mean value, a finite mean energy density and 
satisfies the mixing condition. The distribution of the random solution u{-,t) in the time 
moments t G M is studied. Then the limit 



is established, where /loo is an equilibrium Gaussian measure. For one-dimensional chain of 
harmonic oscillators, this result has been proved by Boldrighini et al. [2] . The convergence to 
equilibrium distribution was proved also for systems described by partial differential equations 
[HI ID], for the crystal coupled to the scalar field [13], and for the Klein-Gordon equation 
coupled to a particle [16] . 

To derive the hydrodynamic equations we apply the special so-called hydrodynamic limit 
procedure in which the notions of hydrodynamic limit and hydrodynamic space-time rescal- 
ings play a central role. Namely, we introduce a small scale parameter e > giving the 
relation between the microscopic and macroscopic space-time scales and consider the family 
of the initial measures {/^O'^ > 0} which satisfies some conditions (see conditions VI and 
V2 in Section 12^2] below) . In particular, we assume that (i) the measures /Iq are locally 



lim jj,t = Id. 
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homogeneous for space translations of order much less than e"'^ and nonhomogeneous for 
translations of order ; (ii) the covariance of /Xq vanishes with distance enough quickly 
and uniformly in e . 

To deduce the Euler equation we use (see [H]) a hyperbolic (or Euler) scaling, i.e., the 
microscopic time and space variables are t = r/e and z = r/e , where the macroscopic time 
and space variables are denoted by r and r , respectively. Given nonzero r e M and r G M'^ , 
we study the distribution /x^^^^^^ of the random solution u{z,t) at the space points close to 
[r/e] G Z'^ and in time moments r/e . Then the limit is established, 

lim/^'/s,^/^ = /i^^^, (1.2) 

where /i^^ is a Gaussian measure (see Theorem 3.3 in JX^ for harmonic crystals in the whole 
space Z'^ and Theorem 2.15 in [17] for harmonic crystals in the half-space ). In particular, 
we derive the explicit formulas for covariance matrix gr,r(^ — z') of the limit measure /i^^ . 
These formulas allow us to conclude that in Fourier transform the matrix function qT,r{(^) , 
9 eT'^ , evolves according to the following equation: 

drf{T,r;e) = iC{e)Vn{e)-Vrf{T,r;e), tGM^ r>0, (1.3) 

where C{9) = ^ q(^q^ ^ ^ ' roughly, Q{9) is the dispersion relation of the 

harmonic crystal (for details, see Section [3.11 below). Here and below dr denotes partial 
differetiation with respect to a time r , Vr/ is the gradient of / with respect to r G M'^ , " ■ " 
stands for the standard Euclidean scalar product in M.'^ (or in ). The equation (11. 3p should 
be considered as the analog of the Euler equation for our model. In |6] , the similar equation 
was deduced in the case d = n = 1 . These results were extended to the harmonic crystals 
in the entire space Z'^ , d > 1 , (see [H]) and in the half-space Z^ = {z E 7^'^ : zi > 0} , see 

The main result of the given paper is the derivation of the equation for the "next ap- 
proximation" to the Euler equation (11. 3p . To obtain the additional term of order e in (11.31) . 
we use a diffusive (or parabolic) scaling, that is we study the distributions of solution u{z,t) 
in time moments of order r/e^. After the appropriate change of variables we derive the 
equation of the Navier-Stokes type 

dJir, r; 6) = tCiO) (viuie) ■ V./(r, r; e) + | tr [V'm ■ V'Jir, r; 9)] ^ . (1.4) 

Here and below V^/ stands for the matrix of second partial derivatives of / with respect to 
r , "tr" stands for trace. The precise statement of the result see in Theorem 13.61 and Corollary 
13. 71 below. In the case d = n = 1 , these results have been obtained in the work of Dobrushin 
et al^. Therefore, in the proof (see Section 6) we use an approach of and tools of pTlfT^ 
developed for harmonic crystals in any dimension. 

In Section [331 we use a scaling of the form t = r/e^ , with k >2, z = r/e, and derive 
the "corrections" of the higher order (i.e. of the order with > 2) to equation (11. 3p . 

Second part of the paper is devoted to the study of the harmonic crystals in the half-space 
Z^ = G Z'^ : 2i > 0} , with zero boundary condition. For such model, we also derive the 
limiting "hydrodynamic" equations (of the Euler and Navier-Stokes types). 
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The paper is organized as follows. In Section 2 we introduce the model, impose the 
main conditions on harmonic potentials and initial measures /Iq and give the examples of 
the potentials and measures /Iq satisfying our conditions. In Section 3 the main results 
are stated. Sections 4 and 5 are devoted to the harmonic crystals in the half-space Z;j[_ . 
Sections 6 and 7 contain the main steps of the proof of results. The technical details of the 
proof are given in Appendices A-C. Appendix D is devoted to locally conserved quantities. 



2 Model I: Harmonic crystals in 

We study the dynamics of the harmonic crystals in Z*^ , d > 1 , 

d{z,t) = - J2 - z'M^'^t)^ zeZ'^, teR, 

z'&i (2.1) 
v{z,{)) =vq{z), v{z,^) =vi{z), zeZ'^. 

Here v{z,t) = {vi{z,t), . . . ,t>„(z,t)) , Vo{z) = (t>oi(z), . . . ,fo„(-2)) G M" , and correspondingly 
for Vi{z) , V{z) is the interaction (or force) matrix, (Vkii^z)) , k,l = 1, . . . ,n . 

Write X(t) = {X%t),X\t)) = iv{;t),v{;t)) and Xq = iXj^,X^) = {vo{-),vii-)) . Then 
(12. ip becomes 

X{t) = AX{t), t e M, X(0) = Xq. (2.2) 
Here A = ( , n ) ' "where V is a convolution operator with the matrix kernel V , 



^ -V ^ 

Vf = X] V{z—z')v{z') . Formally, (12.21) is a linear Hamiltonian system with the Hamiltonian 
functional 

^W = ^El^i(^)l' + ^ 5^ t;oW-n^-^')^o(^'), X = {vo,v,), (2.3) 

where the kinetic energy is given by the first term, and the potential energy by the second 
term. 

Assume that the initial data Xq for (12. 2p belong to the phase space "Hq, , a G M , defined 
below. 

Definition 2.1 T-ta is the Hilbert space of M" x M."' -valued functions o/ 2 G Z*^ endowed 
with the norm ||X||^ = J^zeZ'^ + kP)" < 



2.1 Conditions on the harmonic potentials 

We impose the following conditions E1-E6 on the matrix V . 

El. There are positive constants C and 7 such that ||V(z)|| < Ce~'^'^' for z G Z"^ , where 
||V(z)|| stands for the matrix norm. 

E2. The matrix V{z) is real and symmetric, i.e., Vik{—z) = Vki{z) G M, k,l = 1, . . . ,n , 

zeZ'^ . 
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Let V{0) be the Fourier transform of V{z) with the convention V{6) = ^ V{z)e'^^'^ , 

e eT'^, where T"^ stands for the d-torus M'^/(27rZ)'^ . 

Conditions El and E2 imply that V{9) is a real-analytic Hermitian matrix-valued func- 
tion of ^ G T'*. 

E3. The matrix V{9) is non-negative definite for every 9 eT"^ . 
Let us define the Hermitian non-negative definite matrix, 

n{e) = {v{9))'^^ > 0. (2.4) 

The matrix Q{9) has the eigenvalues < u!i{9) < uj2{9) < . . . < Us{9) , s < n, and the 
corresponding spectral projections Ila{9) with multiplicity To- = trno-(^). The following 
lemma holds. 

Lemma 2.2 (see 07]. Lemma 2.2]). Let conditions El and E2 hold. Then there exists a 
closed subset C T*^ such that the following assertions hold, 
(i) The Lebesgue measure of is zero. 

(a) The eigenvalue 00^(9) , a = 1, . . . , s , has constant multiplicity in T'^ \ . 

(Hi) The following spectral decomposition holds: Q{9) = J2a=i'^o-(^)^^(^) > ^ ^ T'^ \ , 

where Ila{9) is a real-analytic function on T*^ \ C* . 

For 6' G T"' \ , denote by V^u)(j{9) the matrix of second partial derivatives. The next 
condition on V is as follows. 

E4. The functions D„{9) := det (V^a;o-(^)) do not vanish identically on T'^\C* , a = 1, . . . , s . 
Let us write 

Co = {0 G T'^ : det V{9) = 0}, C = G \ : D^{9) =0}, a = 1, . . . , s. (2.5) 

Then the Lebesgue measure of vanishes, cr = 0, l,...,s (see [Hi Lemma 2.3]). Usually, 
the dispersion relations Uk{9) satisfying the condition Uk{0) = are called acoustic. 

E5. For each a ^ o' , the identities C(;o-(6') ± uj^'if)') = const-t for 6' G T'^ \ , do not hold 
with const -t 7^ . 

This condition holds trivially for n = 1 . 

E6. \\V-\9)\\ G L\T'^). 

If Co = , then ||V^^^(6')|| is bounded, and E6 evidently holds. 

Remark 2.3 (i) Instead of condition El we may assume that |V^(-z)| < C(l + 1^1)"''^ with 
an > . However, in this case, we should assume in addition that there exists a set 
JCCT^ such that mes (T'^ \ /C) = and for G /C , G C^JC) , H, G C(/C) , u^{9) ^ , 
det {V^uj^{9)) , Vuj„{9) ^0, a = l,...,s. Note that if condition El holds, /C = \ C 
with C defined in (16. 3p . 
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(ii) Conditions E1-E6 are satisfied, in particular, in the case of the nearest neighbor 
crystal (see [11]) in which the interaction matrix V{z) = [Vkilz))"^ i^-^ is of the form 

{-7fc for \z\ = l, 

2c/7fc + ml for 2; = 0, k = l,...,n, 
for 1^1 > 2, 

with 7fc > and > . In this case, the Hamiltonian functional has a form 

zezd k=i j=i 
Cj = {Sij, . . . , Sdj) , and equation (12.11) becomes 

Vk{z, t) = {-fkAL - ml)vk{z, t), k = l,...,n, 
where Al stands for the discrete Laplace operator on the lattice Z'^ , 

d 

Alv{z) := ^(t;(z + e,) - 2v{z) + v{z - e,-))- 
i=i 

Therefore, the eigenvalues of Q{6) are 



6Ufc(e) = ^27fc(l-cos0i) + ... + 27fc(l-cos^d)+m2, k = l,...,n. (2.6) 

These eigenvalues still have to be labelled according to magnitude and degeneracy as in 
Lemma 12.21 Clearly, conditions E1-E5 hold, and C,, = . If > for any k , then the 
set Co is empty and condition E6 holds automatically. Otherwise, if = for some k , 
then Co = {0} . In this case, E6 is equivalent to the condition uj'^'^{9) G L^(T'^) . Therefore, 
conditions E1-E6 hold if either (i) d > 3 or (ii) d = 1,2 and > for any k . 

Lemma 2.4 (see [77]. Proposition 2.5]) Let conditions El and E2 hold, and a G M. Then 
(i) for any Xq G 1-ia , there exists a unique solution X{t) G C(M, Hq) to the Cauchy problem 



(ii) for any t G M, the operator U{t) : Xq t— )■ X{t) is continuous on Ha ■ 

The proof of Lemma [23] is based on the following formula for the solution X{t) of problem 

X{t)= J2Gtiz-z')Xoiz'), (2.7) 

where the function Qt{z) has the Fourier representation 

g^{z) := F-_]^[exp {A{e)t)] = (27r)-'^ J e''''' exp {A{e)t) dO (2.8) 

Yd 

with 
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2.2 The family of initial measures 

Let £ > be a small scale parameter, {fil,s > 0} be a family of initial measures. To 
formulate the main conditions VI and V2 on the covariance of /Zq , let us introduce the 
complex 2n x 2n matrix-valued function Ro(r, 2;) = (Rq {r, z))} , r G M*^ , 2; G Z*^ , with 
the following properties. 

11. For every fixed r G M'' and i,j = 0, 1 , the bound holds, 

\TCo'ir,z)\<C{l + \z\)-\ zeZ\ (2.10) 

where C is some positive constant, j > d . 

12. For every fixed r G , Ro(r, 6') satisfies 

R°°(r, 6) > 0, Rl\r, 9) > 0, R',\r, 6) = R^°(r, 9)*, 9 G T^. 

13. For every fixed r G M'' and 9 &T'^ , the matrix Ro(r, 9) is non-negative definite. 

gai + ...+aa 



14. For every 9 eT'^ , Rq (■,6') are functions, and the function 



r —7- sup max max 

0^fd «J=0,1 afc=0,l, k=l,...,d 

is bounded uniformly on bounded sets. 



(9r"^ . . . dr'^ 



To derive the equation of the Navier-Stokes type we need an additional condition 14'. 
14'. Ro(-, 9) G Li(M°') , \/9 eJ'^ , and there exist constants C >0 and A^>ci + 3 such that 

sup |Ro(s,^)| < C(l + sgM^. (2.11) 

Here by Ro(s,0) we denote the Fourier transform of Ro{r,9) with respect to r : 

Ro(s, 9) = F,^,[Ro(r, 9)] = j e^'-''Ro(r, 9) dr, s G 9e T'^. (2.12) 

This condition could be weakened (see condition D' in [7J for the case d = n = 1). Instead 
of 14' we may assume that for each 6' G T"' the function Ro(-, 9) admits the representation 



Ro{r,9) = (271)-^ / e-^^X^,ds), (2.13) 



1 

where fi{9, ds) is some Borel (complex-valued) measure on , depending on the parameter 
9 E T'^ . Let |yu|(^^, ■), 6* G T*^ , be the total variation of the measure fi{9,-). This means 
(see [2B, §29]) that for any measurable set E cM'^ , \fi\{9,E) = sup | / f{s)fi{9,ds)\ , where 

E 

supremum is taken over all measurable functions / such that |/| < 1 . We assume that there 
exist constants C > and S > d"^ /2 + 2d + 1 such that for any n G N , 

sup\fi\{9,Kn) <C{l + n)-^, where = {s G M'* : n < |s| < n+ 1}. (2.14) 



Denote by Eq expectation with respect to the measure /ig , and by Q^^{z, z') = Kl{X\z)® 
X^{z')) , z, z' E , z, j = 0, 1 , the correlation functions of /Iq . Assume that Kq{X{z)) = 
and functions Ql^{z,z') satisfy the following conditions VI and V2. 

VI. For any e > 0, z,z' eZ'^, 

\Q'J{z,z') -Ri\ez,z - z')\ < C min \{1 + \z - z' e\z - z'\\, (2.15) 
with the constants C , 7 as in fl2.10p . 

V2. For any £ > and all z, z' E Z'^ , i,j = 0,1, \Q'J{z, z')\ < C{1 + \z - z'\)-^ with the 
constants C and 7 as in f l2.10p . 

Remarks (i) Condition VI can be formulated in the another form (see [HllTl]). Namely, for 
any e > there exists an even integer A^^ such that 

a) for all M e M'^ and z, z' E Im , 

\Q'J{z,z')-Ri^{eM,z- z')\ < Cmin[(l + \z - z'\)-\eNe], 

where C , 7 are the constants from fl2.10p . and Im is the cube centered at the point M 
with edge length N, , Im = {z = (^i, . . . , Zd) E Z'^ : \zj - M,| < NJ2, M = (Mi, . . . , M^)}; 

b) A'e ~ e"^ as £ , with some /3 E (1/2, 1) . 

The formulation of VI in the form fl2.15p is more natural and convenient for our proof, 
(ii) By conditions VI and V2, ^ e'^^-PQ^ {[r /e + p/2], [r/e - p/2]) Iti^{r, 9) as £ ^ 



uniformly in r G M"^ and 9 eT'^ . Here and below [x\ = {[xi], . . . , [xd]) for x G M'^ and [xj] 
stands for the integer part of Xj G , i = 1, . . . ,d . 

2.3 Example of initial measures /ig 

We construct Gaussian initial measures /Xg satisfying conditions VI and V2. At first, we 
introduce matrix- valued functions q^iz) , z eZ'^ , such that q^iz) = for i ^ j , and 

q^{9) = F^^e[q'S{z)]EL\r), &)>0, z = 0, 1. 

Next, we set 

Ri/(r, z) = T{r)qi\z), r G zE (2.16) 

where T G C'^{M.'^) , T{r) > 0, sup sup \D^T{r)\ < C < 00 , |F,^,[r(r)]| < C{l + \s\)-^ 

reR'' \a\<d 

with an N > d . Finally, we put 

Qnz,z') = ^T{sz)T{ez')q^\z - z'), z,z' eZ\ z,j = 0,1. (2.17) 

By the Minlos theorem, for any £ > , there exists a Borel Gaussian measure /ig on "Hq, , 
a < —d/2 , with the correlation functions Ql^{z, z') , because 

E^(ll^ll^) = J2(^ + \z\'rtT[Qf{z,z)+Ql\z,z)] 



= 5^(l + |;.rrTMtr[g°°(0) + go"(0)] 

< C{a, d) tr j (g°°(^) + ql\e)) dO < Ci < oo. 

Let us assume that there exist constants C > and 7 > c? such that 

\qi'{z)\<Cil + \z\)-\ zeZ''. (2.18) 
Then Ql^{z,z') satisfy the conditions VI and V2. 

Definition 2.5 Formally, Gibbs measure g is g{dX) = ^e~^'^'^^'^'^ Y\^^j^ddX[z), where 
H{X) is defined in Ii2.3\) . (5 = , T > is the corresponding absolute temperature. 
We define the Gibbs measure g on Ha, o. < ~d/2, as the Gaussian measure with the 
correlation matrices defined by their Fourier transform as (f^{6) = TV~^{9) , q^^{0) = TI , 
q^^{6) = q^^{d) = , where I stands for the unit matrix in M" . 

Let q^°{e) = V-^{e) , ql\e) = / , e eT'^. Then the functions Rd''(r,2) defined in (^m 
are correlation matrices of the Gibbs measures gr , r G M'^ , with /3 = 1/T(r) . If we assume, 
in addition, that Co = , i.e. det ¥{6) ^ 0, G T*^ , then the bound fl218|l holds, and 
Ql^[z,z') defined in fl2.17p satisfy the conditions VI and V2. 



3 Main results 

Definition 3.1 (i) fif is a Borel probability measure on Tia which gives the distribution of 
the random solution X{t) , fJ^fiB) = fiQ{U{~t)B) , where B G B{'Ha) o-nd t G M . 

(ii) The correlation functions of the measure /i^ are defined by 

Ql^{z,z') =W^{X\z)®X^{z')) =W^{X\z,t)®X^{z\t)), i,j = 0,l, z,z' eZ\ 

where Ef stands for expectation with respect to the measure fif , and X'^{z,t) are the com- 
ponents of the random solution X{t) = {X^{-,t),X^{-,t)) to problem 1^2. 

(Hi) Let Th , h eU^ , be the group of space translations: ThX{z) = X{z — h) , z eTL!^ . For 
r 7^ , r G M*^ , and k > , the measures /i^/^K ^./^ are defined by the rule 

f^T/e'^,r/eiB) = fl';/,.{T[r/s]B) , whcrC B G BiHa). 

Remarks (i) In addition to conditions VI and V2, let us assume that the measures yUg 
satisfy the mixing condition. To formulate this condition, denote by a (A) , A G 'Z'^ , the 
(T -algebra on 'H^ generated by Xq{z) with z E A. Define the Ibragimov mixing coefficients 
of the probability measure fi^ on Ha by the rule 

dist(A B)>r > 
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A measure fi^ is said to satisfy the strong uniform Ibragimov mixing condition if (psi^) ~^ 
as r — 7- oo . Moreover, we assume that Vr G M.'^ , sup^^Q (/^^(r) < C(l + r)^^'^, with the 
constant 7 as in fl2.10p . Note that the last bound on v'e(r) imphes condition V2. 

Then for r 7^ , r G M*^ , the measures /i^^^ converge weakly to a limit measure /x^^ 
on the space Ha , a < — (i/2 . By definition, this means that 

hm I /(X)/i^/,,,/,(rfX) = J /(X)/x^,(c/X) 

for any bounded continuous functional / on Tia ■ Moreover, the limit measure /i^^ is a 
Gaussian measure on Tia (see Theorem 3.3 in ^14j). 

(ii) Let K < 1 . Then for r 7^ , r G M'^ , the following limit holds, /i^/^^ ^ /i^ as e — )■ 
in the sense of weak convergence on Ha, 01 < —d/2. Moreover, the limit measure fir is 
Gaussian, its correlation matrix does not depend on r and has a form (in Forier transform) 

- J2 ^AO) (Ro(r, e) + a(^)Ro(r, 9)C:{9)')n,i9), G T'^ \ C, 

(7=1 

where Co-{0) is defined in fl3.4p . 



3.1 Equation of Euler type 

In this subsection we put k = 1. Let us introduce the matrix qT,r{z) , z E Tj'^ , r G M. 
r G M*^ . In Fourier space, 

^r,M = T$^n.(^^)[$^(/±^a(^^))Ro(r±Va;.(0)r,^^)(/TzC:(^^))jn^(^^) (3.1) 



(7=1 



J2 n.(^^)(M^(r, r;9)+iM'L{T,r;9))U^ie), 9eT''\ 



o-=l 



where Ila{9) is the spectral projection introduced in Lemma [2.21 (iii), 
M;(r,r;^) = '^{Rl{T,r; 9) + C,{9)Rl{T,r; 9)C:{9)), 



M'L{T,r;9) = ^^{C^{9)R'L{r,r;9) - R'L{T,r;9)C:{9)), 
C.{9)= ' 



uj-\9) 



a= l,...,s. 



i?±(r,r;i 



Ro(r + Vu^{9)r, 9) ± Ro(r - Vuj„{9)t, i 



(3.2) 



(3.3) 

(3.4) 
(3.5) 



Theorem 3.2 (^see JT^l Theorem 4-^]) Let the conditions VI, V2 and E1-E6 /ioM. Then 
for any r G M'^ , z^z' EJf' , r 7^ , i/ie correlation functions of measures fi^^^^ converge 
to a limit, 

hmQ^^r/eiir/e] + z, [r/e] + z') = q^^riz - z'). 
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Corollary 3.3 Write f„{T,r\6) = Il^{6)qr,r{G)Ila{0) , cr = l,...,s. Then the function 
f„{T,r;6) satisfies the "hydrodynamic" Euler type equation: 

drUT,r;e) = tC^ie)\/u^ie)-\/rUr,r;e), rGM^ r > 0, (3.6) 

Ur,r;e)\r=o = ^n,{9) (Roir,9) + C^i9)Roir,9)C:i9)^ 11^(9). (3.7) 

Remarks 3.4 (i) Note that qT,r{9) = {ql\.{9))j satisfies the equilibrium condition, i.e., 
ql\i9) = fi2(0)gOO(^), gOi^(^) = -qll{9). Moreover, qi\{9)* = qi\,{9) > 0, g^U^)* = 

(ii) In the case when k G [1,2) , the correlation matrices of measures converge as 

e — i- to the same matrices qT,r{z — z') as in Theorem 13.21 This result can be proved by the 
similar way as Theorem 13. 2[ 

Theorem 13.21 and Corollary 13.31 can be rewritten in the terms of the Wigner matrices. 
Write {V'^v){x) = FiJ^j:[V'^{9)v{9)] , /c G M, and introduce the complex-valued field 

a(x) = (v^/So(x) + iV-'^%iix)^ G C" , X G (3.8) 

with complex conjugate field a{x)* = (l/-\/2) (V^/''fo(x) — iV~^^^Vi{x)) . Define the scaled 
n X n Wigner matrix as 

W%T, r;9)=J2 K/eic^*i[r/e + y/2]) ® a{[r/e - y/2])) . (3.9) 

By properties of /Iq , the following limit exists 

limW^^(0,r;^) = Un'/'R'',\r,9)n'/'' + n-'/^Rl\r,9)n-'/' 

+^fi^/2R°^(r,^)^]-^/2_^^-l/2p^lO(^^^)^l/2^ ^w{0,r;9). (3.10) 



Theorem 3.5 (see fT^ Theorem 3.2]) Let the conditions VI, V2 and E1-E6 hold. Then 
for any r G and r 7^ the following limit exists 

lim Vr^(r, r; ^) = W^{T,r]9) (in the sense of distributions), 

where 

s 

iyP(r, r-9) = Y, n.(^)P^(0, r - tVuj„{9)- 9)T{,{9) . (3.11) 

(7=1 

Remarks (i) The limit correlation matrices qT,r{9) = (^r'!r(^))ij=o (^^e (13. 2p ) are expressed 

by WP{T,r;9) as 

mC((^) = m-'aO) = UW^iT,r;9) + WPiT,r;-9)*), 

C[i]r{d) = -ql'r{9) = -^^{WP{T,T-9)-WP{T,r--9Y). 

(ii) The matrix /o-(r, r; 9) = IIo-{9)Wp{t, r; 9)Ilo-{9) , a = 1, . . . , s , satisfies the energy trans- 
port equation (see [m p. 656]) 

drfAr,r;9) + Vu,{9)-VrfAr,r;9) = 0, r > 0, r G 

with the initial condition f^{T,r;9)\r=o = Il^i9)W{0,r;9)Il^{9) , r G M*^ . 
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3.2 Navier— Stokes equation 

In this subsection we study the behaviour (as £ — )■ ) of the correlation functions of /i^/^^; 
for K = 2, and obtain the next term of the decomposition in e to equation fl3.6p . Let us 
introduce the matrix , ^ G Z'^ , r G M'^ , r 7^ , e > , which has the following form 

(in the Fourier transform) 

1 

^rM = 7 E n-(^) f ± ^a(^))Aj,(r, r; 0)(/ T ^C7:(^^))j n.(0), G T'^ \ (3.12) 



4 

where matrices CaiO) are defined in (13.40 . 

A%^{T,r-e) = jRoir±Vu^{e)T/6-x,e)K^{T,x,e)dx, (3.13) 

ir±(r,a;,^) := F-^Je^F'^"/')^-^^'"'^^'))^], a; G R'^ (3.14) 
(see also formula (16.251) ). 

Theorem 3.6 Let conditions 11-14', VI, V2 and E1-E6 hold. Then for any r 7^ , r G 

M.'^ , z,z' G 'Z'^ , the correlation functions of measures /i^/£2 ^./^ have the following asymptotics 

lim (Q,,,/,2([r/£] + z, [r/e] + z') - q%{z - z')) = 0, (3.15) 



^-1 \f.e 



where the matrix q^^ri^) = -^g^zl^r.rl^)] ^■^ defined by ^3.12\} . 

We omit the proof of this theorem since it can be proved by the similar technique as 
Theorem 15. 6^ below. 

Remark. Note that (f^^{0) satisfies the equilibrium condition (see Remarks 13.41 (i)). 
Set T = et. It follows from formulas f lXT^ -f lXTi]) that 

A%^{et, r- e) |,=o = Ro(r ± Vuj,{e% 6). 

Hence, g|j^(6')|e=o = Qt,r{0) , where qt,r{6) is defined in (13. ip . Denote 

ll,...,lfc = l K L K 



Corollary 3.7 Let r G M'^ , t G M, Q ^T'^\C^ . It follows from formulas [SJE )- ^3l4\ ) that 



for each a = l,...,s, the matrix-valued function f^{t,r;9) = Il^{9)q^fj.{9)Ila-{6) evolves 
according to the following (Navier-Stokes type) equation 

dj',{t,r;e) = tC^{e)(Vu^{e) ■ V, + | W(^) ■ V^)/^(t,r;^), teR, rem", (3.17) 



with the initial condition l\3. T] ). 
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3.3 Corrections of the higher order 

To obtain the "corrections" of order , e'^ ,. . . to equation fl3.17p . it is neccesary to study 
the behaviour (as e — )■ ) of the correlation functions of measures //^ ^.^^ with n> 2 . 

Theorem Z.S Let k>2 , r ^W^ , t > Q . Then V^, z' elJ^ , 

Qe,r/e^(\r/e] + z, [t / s] + z') - ^^^(^ " ^ 0, 0. 

The matrix q^'^iz) (k = 2,3, . . . ) has the following form (in the Fourier transform) 

^r'M = -Y^ ii.{e) [Y,{i ± ^cM)At^{r, r; e){i T ^c:{e))\ u^{e), 

where 



Kf;^{T,x,e) = F;X. 
Here, by definition, 



2! A;! ' 



x,y e R'^. 



dOi, . . . dOi, 
«i,...,jfc=i 



Note that the matrix (?^'^(^) coincides with q^^i^) from formula f l3.12p . The proof of 
Theorem 13.81 are similar to the proof of Theorem 13.61 

Set T = e'^-H, k>2. Then dtAf'^{e^~H,r; 9) = ±P^{9,dr)Af^{e''-H,r;9), where 

P^{9,dr) := ^^VPa;.(^) ■ 

(see notation (13.161) ). Therefore, the matrix Ilu{9)q'^^iLi^ ^{9)Ilcr{9) (denote its by f^(t,r; 9) ) 
is the solution of the following equation 

dtft{t,r-9)=iC„{9)P^{9,d,)ft{t,r-9), t > 0, rG< (3.18) 

with the initial condition (13. 7p . 

Denote by f^^ , i,j = 0,l, the elements of the matrix f^(t, r; 9) . Then fl^ = uj1{9)f^^ , 
fl^ = — , and the equations (13.180 can be rewritten in the form 

= {-^M9))P,'i9, 9.)/°\ = tu^{9)P^{9, 9.)/^- 
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4 Model II: Harmonic crystals in the half-space Z'J. 

We study the dynamics of the harmonic crystals in Z'^ , d > 1 , 

M(z,t) = - ^ zeZ^, t G M, (4.1) 

with zero boundary condition, 

uiz,t)U=o = 0, (4.2) 

and with the initial data 

u{z,0) =uo{z), u{z,0) =ui{z), zeZ%. (4.3) 

Here Z'J_ = {2; G Z*^ : > 0} , z = {—zi,Z2,---,Zd)- For convenience, we assume that 
Uo{z) = Ui{z) = for zi = . 

Write F(t) = iY%t),Y\t)) = («(-, t), «(-, t)) and Yo = (Y^^X) = (%(■), • Then 
fl4.ip - fl4.3l) becomes the evolution equation 

Y{t) = A+Y{t), t G M, z G Z^ Y\t)U=o = 0, Y{0) = Yo. (4.4) 
Here A+ = ( ° I] with V+u(z) := ^ (\/(z - z') - V{z - z'))u{z') . 



V+ 
Let us assume that 

V{z) = V{z), where z = (-^1,^), z = {z^, . . . , z^) e (4.5) 

Then the solution to problem (14. 4p can be represented as the restriction of the solution to the 
Cauchy problem (12. ip with odd initial data on the half-space. More exactly, assume that the 
initial data Xo{z) form an odd function with respect to zi G Z\ i.e., let Xo{z) = —Xo{z) . 
Then the solution v{z,t) of (12.11) is also an odd function with respect to zi G Z^ . Restrict 
the solution v{z,t) to the domain Z'^ and set u{z,t) = v{z,t)\z-^>o ■ Then u{z,t) is the 
solution to problem (14.10 with the initial data Yo{z) = Xo{z)\z-^>o . 

Assume that the initial data Yo for (14. 4p belong to the phase space ^^0;,+ , a G M , defined 
below. 

Definition 4.1 ?^q,,+ is the Hilhert space of x R" -valued functions of z ^ Z"^ endowed 
with the norm \\Y\\l^^^ = J2zez^ \Y{z)\'^{l + I^H" < 00. 

In addition, it is assumed that the initial data vanish ( Fo = ) at = . 

Lemma 4.2 (see fTSi Corollary 2.4]) Let conditions El and E2 hold. Choose some a G R . 
Then (i) for any Yo G 'Ha.+ , there exists a unique solution Y{t) G C(R, 7^0,+) to the mixed 
problem (^T^; 

(ii) the operator U+{t) : Yq Y(t) is continuous on 710,+ ■ 

Lemma 142) follows from Lemma since the solution Y{t) of ( 14. 4p admits the representation 

Y{t) = GtA^^ z')Yo{z'), zeZ%, (4.6) 

where Qt^+^z, z') := Qtiz - z') - Gtiz - z'), (4.7) 
and Gt{z) is defined in (12.80 . 
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4.1 The family of the initial measures 

Let us introduce the complex 2n x 2n matrix-valued function i?(r, x, y) = {R^^{r, x, y))\j^Q , 
r e M*^ , X, 7/ e , with the following properties (a)-(d). 

(a) R{r,x,y) = for xi = or ?/i = . The n x n matrix- valued functions R^^{r,x,y) 
have the form 

R'^{r,x,y) = R'^{r,xi,yi,x -y), where x = (xi, x), y = {y^,y), z,j = 0,l. 
Moreover, uniformly in r G M'^ , 

lim R''{r,y, + zi,yuz) = B!^{r,z), z={zuz)eZ\ z,j=0,l, (4.8) 

where the matrix Ro(r, 2;) = (Rq '^))ij=o satisfies conditions 11-14 (see Section [2^ . 

(b) For every fixed r G M°' and i, j = 0, 1 , the bound holds, 

\R^{r,x,y)\<C{l + \x-y\)-\ x,yeZ'i, (4.9) 
with the same constants C and 7 as in (12.101) . 

(c) For every fixed r G M'^ , the matrix- valued function R{r, x, y) satisfies 

i?"(r,-,-)>0, R'^{r,x,y) = {R^%r,y,x)f, x,yGZ^. 

(d) For every x,y E Z'^ , R'^{-,x,y) , i,j = 0, 1 , are functions. 

To derive the equation of the Navier-Stokes type we need the additional condition (d') 
in the case when d > 1 . 
(d') Let d > 1 . For every fixed x,y E Z'^ , 

sup \R{ri,f,x,y) — R{0,f,x,y)\—>0 as ri — )■ 0. (4-10) 

f=(r2,...,rd)eM''-i 

Moreover, we assume that for every fixed x, ?/ G Z'^ , 

|^(0,s,x,i/)| < C(l + |s|)-^, sgM'^-\ (4.11) 

with N > d + 2 . Here by -R(0, s, x, y) we denote the Fourier transform of -R(0, f, x, y) w.r.t. 

f G R'^-^ , 

^(0,s,x,?/)= y e^'"^/2(0, f, X, y) rff, seR'^-\ 

Let {/Uq'^ > 0} be a family of initial measures on 'Ha,+ and Eq stand for expectation 
with respect to the measure /Iq . Assume that E,q{Y{x)) = and define the covariance 
Q'Jix.x') = Eg(F^(x) (g) Y^{x')) , x,x' G , z,j = 0, 1 . 

The family of measures {/Ug,e > 0} satisfies the following conditions VI' and V2'. 
VI'. For any e > 0, x, x' G Z^J. , 

\Q'J{x,x') - R'^{ex,x,x')\ < Cmin[(l + |x - x'])"^, e|x - (4.12) 

with the constants C and 7 as in (14.91) . 

V2'. For any e>0 and all x,x'gZ^, i,j = 0,1, |Q*^(x,x')| < C{1 + \x - x'\)-^ with the 
constants C , 7 as in (14.91) . 
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5 Main results in the half-space 

Definition 5.1 (i) fi^ is a Borel probability measure on Ha,+ which gives the distribution 
of Y{t) , fiKB) = ij,l{U+{-t)B) , where B e B{na,+) and teR. 
(a) The correlation functions of the measure /xf are defined by 



Here Y^{x^t) are the components of the random solution Y{t) = (Y^{-,t),Y^{-,t)) to the 
problem ^■4\)- 



5.1 Euler limit 

At first, we introduce the matrix gr,r{z) , -2 G Z'^ , r G M*^ , r 7^ , by tlie Fourier transform, 

9rM = T E [ ± ^cM)Mr ± va;.(0)r, e)x%, {e){i T ^c:{e))\ n.(^), 

(T=l ± 

(cf dSH)), where G T"^ \ , 

X%,{0) = (1 + sign(ri ± d,ujMr))l2. (5.1) 



Theorem 5.2 (see Theorem 2.10 in JT^/; Let conditions (a)-(d), VI', V2' and E1-E6 
hold. Then for any r 7^ , r G M'^ with ri > , the correlation functions of /i^/^ ^./^ 
converge to a limit, 

limQ,,,/,([r/£] + z, \r/e] + z') = Q^A^, z'). (5.2) 
Here z, z' eU^ if ri > , z, z' e Z'l if n = , 

^ \ 9tAz-z') - grAz- Z') - grAz- Z') + grAz- Z'), Ti = 0, 

where r := (— ri,f) , f = (r2, . . . ,rd) , and the matrix gr,r{z) is defined above. 

Corollary 5.3 Let r e R'l = {r E : ri > 0} and r > . For each a = 1, . . . , s , the 
a -band of cir,r{0) , 6* G T'^\C* , satisfies the following "hydrodynamic" (Euler type) equation: 

drfAr, r; 6) = tCAO)VioAO) ■ V./.(r, r; ^), r G ^ > 0' 

with the initial condition ( [j*. ?| j (if t = ) and with the boundary condition (if ri = ) 
expressed by Rq . In particular, if Ilo{r,6) = Ro(r, ^^) , the boundary condition has a form 

M,,=o = ^A0)1{PI + cAe)P+c:{e) + iCAe)Pi - iPic:{e))iiAe)^ f e r''-\ r > o, 

where stands for the 2n x 2n matrix-valued function, 

Pl = \ (Ro {\diUjAO)\r, f + V-eUjAO)r, 6) ± Rq (|9ia;<,(^) |r, f - Ve-a;.(0)r, 9) ) . 
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Introduce the complex- valued field (cf. (13. 8p ) 

a+(x) = -^((Vy'uo)(a:) +^(V;'/V)(a;)) G , a; e Z^ (5.4) 

where 

(V^n)(x) := J2 (^^(^ - ^) - ^''(^ - with l^'=(2) := F^I^^ (^V^'=(0)' 

Let us introduce the scaled n x n Wigner matrix (cf. (13. 9p ) 

W^{T,r;e)=Y,e''-'K/e{<{[r/e + ym)®ci^{[r/e-y/2])^, r G 

where a+(x) is given in (15.40 . By conditions VI' and V2', for fixed r G M!^ , 

limiy"(0,r;e) = W(0,r;e), 

£-5-0 

uniformly on G T'^ \ C* , where W{0, r; 6) is defined in (I3.10p . 
We also define the limit Wigner matrix as follows (cf (13. lip ) 

W^{T,r;6) = I '^7^ (5.5) 

I j2^Ad)wio,-n+Tdiu^ie),f-TVou^ie);e)u^ie),ii n<Tditu^ie), 

V a=l 

where 6 = {-6^, 9) , 9 = {02, ...,9,), VeU^{e) = {d2Uj^{e), d,uj^{e)) , d, = 8/89^ . 
Theorem 5.4 (see Theorem 2.12 in fl^) Let conditions VI', V2' and E1-E6 hold. Then 



for any r G M!^ and r > , the following limit exists in the sense of distributions, 

\imW^{T,r;e) = Wl{T,r;e). (5.6) 

Corollary 5.5 (see Corollary 2.13 in fF^ ) Denote by W^{T,r;9) , a = 1, . . . , s , the a -th 
band of the Wigner function W^{t, r; 9) . Then, for any fixed a = 1, . . . ,s , is a solution 
of the "energy transport" equation 

drWP{T,r;e) + Vu^ie)-VrW^{r,r;e)=0, r > 0, r G 

with the initial and boundary conditions 

WP{T,r;e)\^=o = W^{0,r;e), r G 

W^iT,r;e)\r,=o = K{r,f-9), f eW'~\ r > 0. (5.7) 

Here Wcj{Q,r\9) is the a -th band of the initial Wigner matrix W{0,r;6) (see Ii3.10\) ). 

W^iO, -rdiiu^ie), f - rVgiO^iey, 6), if diu^ie) < o, 
W^{0,Tdiu„{9),f-rVeUJ„{9y,9), if diu^{d) > 



In particular, if we assume that Ro(r, 6) = Ro(r, 6) for r eR"^, 6 eT'^ , then the boundary 
condition ( [5. 7| ) can be rewritten in the form 

ws{T,r;e)\r,=o = w^{o,T\dioo^{e)\,f - TVgoo^{ey,e), f g m^-\ r > o. 
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Remark (see [T71 Theorem 2.15]) Let the measures /Xq satisfy the mixing condition (of the 
Rosenblatt or Ibragimov type). Then for r 7^ , r G M"^ with ri > , in the sense of 
weak convergence on 1-La,+ , liin yU^/^ = /^^^ ■ The measure /x^^ is a Gaussian measure 

on 'Ha,-!r 1 which is invariant under the time translation U^{t) . fi^^ has mean zero and 
covariance Qr,r{z,z') defined by (15. 3p . 



5.2 Second approximation 

In this subsection we treat the main result of this paper. Let us introduce the matrix g^ri^) , 
z & 'Z'^ , rGM'^, T ^ , £:>0, which has the following form in the Fourier transform (cf 
dffl)) 

9r,M = \i2 n-(^) [ ± iC^{9)) Afjr, r- 9){I t ^C:{e))] U^{9), Oe^X C,, (5.8) 
where the matrices C„{9) are defined in (13.41) . 

with the matrix-valued functions Af^{T,r;6) from (I3.13P and xtn from (15. ip . 

Theorem 5.6 Let conditions (a)-(d'), VI', V2' and E1-E6 hold. Then for any r 7^ 
and r G M"' with ri > , the correlation functions of measures /Xr/e^^r/e have the following 
asymptotics 

hm (Q,,./,2([r/£] + z, [r/e] + z') - Q%{z, z')) = 0, (5.9) 
where z, z' ^if' if ri > , z, z' e Z'^ if n = , 

y'\ = \ ^^''■'^^ " ^'^ ^ ^'-^^^ ~ + ^'-^^^ ~ ^1 > 

^^'^^ ' ' I gU^ - z') - g%{z - z') - g%^X^ - z') + ^7^,,(5 - F), z/ n = 0, 

g^^^^iz) = F.^J^g^^^M] and g%{e) defined by (EJ). 

Remark. The matrices qrr(^) and cf^^{6) satisfy the equilibrium condition (see Remarks 

Ea](i)). 

Set T = et . In this case, A^^(et, r; 6')|e=o = Ro(^ i Vu;a-(6')t, ^)xfri(^) 5 ^"^^i then, 
9lt,r{^)\e=o = gt,riO). Therefore, for ri > , ci%^^.{9)\e=o = (lt,r{0) , where qtA^) = 
F^^e[qt,r{z)] and qt,r{z) from (jSJD- 

Corollary 5.7 Let ri > anc? t > 0. T/ien for each a = l,...,s, the matrix-valued 
function = F^{t, r; 9) = Il(j{9)ql^ ^{9)Il„{9) , 9 G T"'\C* , evolves according to the following 
mixing problem 

dtF„ = iCA9) [vuj,{9) ■ VrF^ + I W(^) ■ VlF}j , n > 0, t > 0, 

F^t=o = ^n.(^)(Ro(r,^) + a(^^)Ro(r,^)c:(^))n,(^), n > 0, 

F.|.,=o = ^^n,(^)(/±2a(^))A;,(£t,r;^)|,,=o(/T^C7:(0))n.(^^), t>0. 



± 
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6 Convergence of correlation functions 

6.1 Bounds for initial covariance 

Definition 6.1 By P = £p(Z°') ® R" (by = £p(Z^) O M") , where p > 1 and n > 
1, denote the space of sequences f{z) = {fi{z),...,fn{z)) endowed with norm \\f\\ep = 

{E.ez^\fi^)\'y\ respectively, = (E.ez^ l/WI")'^' ■ 

The following lemma follows from condition V2. 
Lemma 6.2 Let condition V2 hold. Then, for i,j = 0, 1 , the following bounds hold: 

J2\QHz,z')\ < C<oo forallzeZ'l, 
\Q'J{z,z')\ < C <oo for all z G Z^. 

Here the constant C does not depend on z, z' G Z'J_ and e > . 
Corollary 6.3 By the Shur lemma, it follows from Lemma \6. 2\ that 

\{Q,{z,z'),<!>{z) (^^{z')) + \ < C\m,2jnil, forany <i>,^ e il, 
where the constant C does not depend on e > . 

6.2 Stationary phase method 

By dal]) and ([2l]) we see that QtiO) is of the form 

A / cosQt sin fit Q'^ \ ,„ 

^*(^)=^-sinfitl] cos^]^ J' ^^-^^ 

where VL = Q{6) is the Hermitian matrix defined by (12.41) . Hence, by Lemma [2.21 (iii), by 
formulas cosw<,(^)t = (e*'^'^* + e-^'^"*)/2, sinw^(^)t = (e^'^'^* - e-*'^"*)/(2i) and by ([314]), the 
matrix Qt{x) can be rewritten in the form 

g^^a;)=J2 I e-'''-'e^'^"^'^'c'^{e)dd, (6.2) 

where c'^iO) := TiCaiO)) /2 . We are going to apply the stationary phase arguments 

to the integral (16.21) which require a smoothness in 6 . Then we have to choose certain smooth 
branches of the functions 0^(6*) and oj^iO) and cut off all singularities. First, introduce the 

critical set as 

C = C„ UC U (C, U {» e T" \ C. : %P = 0} U {« e C. : = O}) . (6.3) 

(T=l i=l « ^ 
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with C* as in Lemma \T2\ and sets Cq and defined by fl2.5p . Obviously, mesC = (see 
[TT| Lemmas 2.2, 2.3]). Secondly, fix an 6 > and choose a finite partition of unity, 

K 

f{e) + g{e) = i, g{9) = Y,9k{e), eer, (6.4) 

k=l 

where f,gk are non- negative functions in C^{T'^) , and 

supp f c{eeT'^: dist(^, C) < 5}, supp gkC{eeT'^: dist(^, C) > 5/2}. (6.5) 
Then we represent Qt (x) in the form Qt (x) = (x) + Qf (x) , where 

Qli^) = FeUfid) QKx) = FiUm (6.6) 

By Lemma 12.21 and the compactness arguments, we can choose the supports of gk so small 
that the eigenvalues Ua{9) and the amplitudes c^(^) are real- analytic functions inside the 
supp gk for every k . (We do not label the functions by the index k to simplify the notation.) 
The Parseval identity, fl6.ip . and condition E6 imply 

\\gI{-)\\% = c J\gt{e)\''\f{e)fde<c J \gt{e)f de ^ o s ^ o, 

Jd dist(0,C)<(5 y-j 

\\gf{-)\\% < c j i^<(^)prf0<Ci<oo, 

Yd 

uniformly in t G M. For the function Qt{x) , the following lemma holds. 

Lemma 6.4 (see fL^ Lemma 4-'^]) Let conditions E1-E4 and E6 hold. Then the following 
hounds hold. 

(t) sup,,^.|^f(x)|<Ct-'^/2. 

(a) For any p > 0, there exist Cp,jg > such that |^f(x)| < Cp{\t\ + |x| + 1) ^ for 

\X\ > 'Jgt . 

6.3 Proof of Theorem 5.6 

The representation (14.61) yields 

for any t G M . It follows from condition (14. 5p and from formulas (12. Sp and (12. 9p that 
Qt{z) = Qt{,z) with z = {—zi,Z2,...,Zd). In this case, by (14. 7p . the covariance Qs^z.z') 
can be decomposed into the sum of fourth terms, 

Qe,t{z, Z) = Se,t{z, z!) - Se,t{z, z') - Se,t{z, z) + Se,t{z, z), Z, Z G Z^, 

where 

Se,t{z,7^) := Qt{z-x)Qe{x,y)Qt{^ -yY . 
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Proposition 6.5 Let r G M°' and z,z' E TL^ . Then 

^Um (^,,./,2([r/£] + z, [r/e] + z') - g%{z - z')) = 0, (6.8) 

where g^^.{z) is defined by ( 15. <^) . 

This proposition implies Theorem 15.61 Indeed, let ri = . Then f/e = r/e, and for 

z, z' eZ'l, 

QsA[r/e] + z, [r/e] + z') = SeAV/e] + z, [r/e] + z') - Se,t{[r/e] + 5, [r/e] + z') 

-S,,{[r/e] + [r/e] + z') + S,,{[r/e] + 5, [r/e] + 5'). 

Therefore, convergence fl6.8p imphes f l5.9p . 

Let ri > . In this case, the matrix-valued functions S^^r/e'^i\f /^] + 5? ['^Z^] + z') and 
'S'e,T/e2([r/£:] + z, [f/e] + 5') vanish as e — )■ +0 , and 

5,,,/,2([f/£] + z, + z') - gl -{i - 5') ^ 0, £ ^ +0. 

It can be proved similarly to Proposition 16.51 
Proof of Proposition 16.51 Let us denote 

n ,A - S Qei^^y) for ^'2/ ^ 
'^^ \ otherwise 

The partition (16. 4p . Corollary 16.31 and the bound (16. 7p yield 

SeAz,z')= G!{z-x)Qs{x,y)g!{z'-yf + o{l), 

where is defined in (16. 6p . o(l) — i- as 5 — i- uniformly in t G M and z,z' G Z'^ . In 
particular, setting t = r/e^ , z = [r/e] + 1 and z' = [r/e]+ p we obtain 

S,^,/,.{[r/e]+l,[r/e]+p) = S',/eA[r/e] + I - x)QAx,y)g'^/^,{[r/e] + p - yf + o{l) 

= E ^.V(/ + x)Q.([r/e]-x,[r/e]-y)^f/^,(p + y)^ + o(l).(6.9) 

Let c = 7g + max(|Z|, \p\) with 7^ from Lemma [6.41 Then Lemma [6.41 (ii) and condition V2' 
imply that the series in (16. 9 p can be taken over x,y E [— cr/e^, cr/e^]'^ . 

By definition, the function R is equal to 



R{r,x,y) 



R{r, X, y) if x,y E Z'^, 
otherwise 



The asymptotics of S^^^i^2 is not changed when we replace Qe{[r / e]—x,[r / e]—y) in the r.h.s. 
of dnSD by R{. . .) = Rie[r/e] - ex, [r/e] - x, [r/e] - y) , i.e., 

Se,r/A[r/e]+l, [r/e]+p) = Yl ^^V^^ + ^)^(- ' -^^r/e^^P + + ^(1)- (6.10) 
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Indeed, by Lemma 16.41 (i) and condition fl4.12p , 

I E + ^) {Qe{[rle] -X, [r/e] -y) - R{. . + V? 

x,y^\—CT I fir I e^Y" 

< C min[(l + \x — y\)~'^, £b\x — y\], 



where the last series is order of . This goes to zero as e — )■ , since 7 > c? . Using 

Lemma [6.41 and properties of i?, we can take the series in f l6.10p over G Z"^. Hence, 

S,^,/,.{[r/e] + l,[r/e]+p)= S'^/A^ + x)R{. . .)g'^^^,{p + yf + o{l), e ^ 0. (6.11) 



Let us spht the function R into the following three matrix functions: 

R^{r,x,y) := ^Ro(r, a; - |/), (6.12) 

R-{r,x,y) := ^Ro(r, x - y) sign(yi), (6.13) 
R^{r,x,y) := R{r,x,y) — R'^{r,x,y) — R~{r,x,y). (6.14) 
Next, introduce the matrices 

S:,r/e-^ = Sl^;A[r/e] + l,[r/e]+p) 

= E ^r/sA +x)R''ie[r/e] - ex, [r/e] - x, [r/e] - y)g',/AP + vf^ (6-15) 

for each a = {+, — , 0} and split Se,T/e^ into three terms, Se,T/e^ = SA/e^ + S~^/^2 + SA/£2 ■ 
The convergence (16.81) results now from the following three Lemmas I6.6H6.8[ since (see (15. 8p ) 
9r,M = (l/2)g;,(2) + il/2)f^^^^.{z) , where g^^,(z) is defined by (IXT^ and f^/z) is defined 
in Lemma [6.71 

Lemma 6.6 lim(S'+^/^2([r/£] + /, [r/e]+p) - (l/2)g^,,(/ - p)) =0, l,peZ'^ , where q'^^^{l) 
is defined by h3.1l!\) . 

Lemma 6.7 lhn(S'-^/^,([r/£] + /, [r/e] + p) - (l/2)/,-^,(/ - p)) = , /,pGZ'^. 

Here f^A) stands for the matrix which is defined similarly to q^r{l) but with the matrix 
Af^^{T, r; 6) sign(ri ± diUJ„{9)T / e) instead of Af^^{T, r; 9) in the r.h.s. of ^3.12\) . 

Lemma 6.8 lim5°^/^2([r/e] + /, [r/ej+p) = , l,peZ'^. 

The proofs of Lemmas 16.71 and 16.81 see in Appendices A and B, resp. 

Proof of Lemma 16. 6[ Step (i): By (16.120 and (16.150 . the function S^i^i can be represented 

as 

^tr/e^ = \Y. ^r/eS + ^) E / s] ' sx^y - x)^f/^,(p + yf. 
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Using the Fourier transform and the Parseval identity we rewrite ,0 as 

Change variables in ( I6.16P : {9, 9') {9,ip) , ip = 9 — 9' and rewrite in the form 

^;./,. = (2^)-'^ / e-'^'-^>'lMO'r/A(^r d9, (6.17) 



d 



where Ie{9) stands for the matrix- valued function, 

IM = J^C^^y 1 e'^'+^^-^Ql,A(^-^)d^Me[r/e]-ex,9). (6.18) 



Step (ii): For simplicity of the proof, let us assume that 'Ko{r,9) satisfies condition (12. lip . 
Under the more weakened condition fl2.14p the proof is given in Appendix C. 
We apply the Poisson summation formula (see, for example, [2^) and obtain 

e^^-^Ro(£[r/£] - ex, 9) = e'^ ^ Ro(-^/£ - l-nnje, 9y^'l'^-^, (6.19) 

where Ro(',^) stands for the Fourier transform of Ro(',^) (see f l2.12p ). We substitute the 
last expression to f l6.18p and obtain 

h{9) = (27r£)-^ f e^^^^VI-\)-^Ql^^, _ ^ Ro(-<^/£ - 27™/^, 9)) + o(l). (6.20) 



Jd 



Condition 14' implies that 



£-'^^Ro(-¥?/e-27rn/£,^) ^ as e ^ (6.21) 

uniformly in ^ € T'^ , G {-tx.txY . Hence (ICTl) and ( lOTl) yield 

/e(^) = (27r£)-'^ y e^'^-('+['-/^l)^f/^,(^-<^)Ro(-<^/£,e)d<^ + o(l). 

Step (Hi): Change variables (p — )■ —eip in the last integral and obtain 

W) = (277)--^ J e-'^-('+['-/^])^f/^,(e + £<^)Ro(<^, 9) + 0(1). 

[— 7r/£,7r/e]'* 

Now we use the following representation for 0^(9) (see (16.21) ): 

^tW = EE^^^""^'^*^^(^)' ^^ere h^i9) = g{9)c^{9)=gi9)^-^^^n^i9), (6.22) 
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and rewrite le{0) in the form 



"'^ l-TT/e,n/e]'' 



since e"*"^ ' - 1 = 0{e) and e"*^-"!''/^] - e"^'^-'' = 0{e) by (|21T|) . 
S'te^} ('zwj: Let us take the Taylor sum representation for IjJo-{0 + eip) : 

+ £(/.)r/£2 = io,{e)T/e^ + if ■ Vuj^{e)T/e + {t/2)^ ■ H^{e)ip + 0{e\ipf), (6.23) 
where H„{9) denotes the matrix V'^uJa{0) . Hence, by (12. lip . 

IM = K{dy^'^^^'^^''"Al^{r, r- 9) + o(l), e ^ 0. (6.24) 

Here (see formulas fl3.13p and (I3.14p ) 

AfjT,r;e) = (27r)-'^y"e-'^-('^^^"'^(^)"/^)e±'("/2)^-^'^(^)^Ro(yp,^)c/v? 



j Ro(r T ^i^Mr/e - x, e)K^{T, x, 6) dx, 



and K^{t, x,6) ( r > , x G M*^ , 6 E T'^\C) stands for the following matrix- valued function 



g±«7rs/4 ^=fi/(2T)x-H~'^{0)x 



K^{t,x,9) = [e±*{-/2)^-/^.W^] = " " , (6.25) 



where s denotes the signature of the matrix H„{9) . By definition, the signature of the non- 
degenerate symmetric matrix A means the signature of the quadratic form with this matrix 
(or the difference between the sums of positive and negative eigenvalues of A). Formula 
fl6.25p follows from the following equality (see, for example, Ramanujan's integrals in [2^ 
section 7.5] or [23, §3]): for any a > , 



hoo 



V a 

Step (v): We substitute flCTp in (leTTI) . apply the decomposition (16:22|) to , ^{0)* and 
obtam 



Str/s^ = (2vr)-'^^y"e-^('-^)-^(g/.^(^^)e±^-(^W^^A,^,(r,r;^ 

(^e±*"-'(^)"/^'/i;^,(e)*)rf5 + o(l), e^O. (6.26) 



Yd 
X 



<t',± 
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Therefore, to find the asymptotics of S^^i^2 it suffices to study the behaviour as £ — )■ of 
the following integrals 

= (27r)-'^i j e-*('-P)-^e'("-(^)±"-'(^))"/"'/i;(0)A-^(r,r;0)/i;^,(0)*rf^, (6.27) 
Jd 

cr, cr' = 1, . . . , s . Note that sup \A~^{t, r]6)\ < C < oo by condition 14'. Hence, the func- 

eeJ<i\c 

tion h~{6)A~^{T,r;6)h^,{d)* G L^{T'^) by condition E6. Therefore, the oscillatory integrals 
with ujfj{6) ±u„i{9) ^ const-t vanish as e — ?■ by the Lebesgue-Riemann theorem. Further- 
more, the identities uJcr{0) ±ijJa-'{0) = constj- in the exponent of fl6.27p with constj. 7^ are 
impossible by condition E5. Hence, only the integrals with 00^^(6) — uj^jiiO) = contribute 
to the integral (I6.27P since 0J^{9) + uj^i{0) = would imply uja{0) = oJa'^O) = which is 
impossible by E4. Therefore, if a a' , lf^i{e) = o(l) as £ ^ . If a = a' , I^„{e) = o(l) , 

4- (e) = (27r)-'^i J e-^'^^-^>'h-{e)A-^{T,r;e)ht{erde + o{l), e 0. 

Jd 

Finally, using (K^ and the equalities h^{9) = g{9)n^{9){I t iC^iO)) /2 and flXT^ . we 
obtain 

= (2^)"'^ / E E ^5 ^)^^ (^)* de + 0(1) (6.28) 



7 Convergence of Wigner matrices 

Here we prove Theorem 15. 4[ Theorem 15.61 implies that for any fixed r G M!J[. , r 7^ , and 
y G (2Z)'^ , the following convergence holds, 

lhnE^/,(a;([r/£] +y/2) ® a+([r/£] -y/2)) = W^(r,r;y), (7.1) 
where in the Fourier space one has 

+zn'/^ql],{9)n~'/^ - zr]-i/2qio (^)^i/2 j _ y^P^^^ ^.Q^^ reWl, (7.2) 

by formulas (15. 5p and (15. 3p . Therefore, convergence (15. 6p follows from (17. ip . (17. 2p and the 
following bound: 

sup sup sup |Ef(a^(z) ® aj^{z'))\ < C < 00. 

e>0 tm z,z'eZ<* 

The proof of this bound follows from Lemma 17.11 
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Lemma 7.1 Let conditions \ 2'' and E1-E3, E6 hold and let a < —d/2 . Then the following 

bound holds: supsup sup < C < cxd. 

£>o tm z,z'&'l_ 

Proof. The representation f l4.6p gives 

y,y'&% fc,'=o,i 
= {Q,{y,y'),^\{y,t)®^i,{y',t)) + , 

where ^\{y,t) is given by 

^'Ay.t) = {gi'4z,y),gil4z,y)^ 

= {gf{z-y)~gf{z~y),gi\z-y)-gi\z~y)), t = 0,l. 
The Parseval identity, formula f l6.ip . and condition E6 imply that 

\m-M = i2n)-'l me,t)fde<c I {\gfm' + \gi\9)\') de<Co<oo, 

where the constant Cq does not depend on z e Z*^ and t G M. Corollary gives now 

mi^,^')\ = \{Qs{y,y'),<^l{y,t) ® Kiy',t))+\ < c\\<i>i{;t)\\ii WKi-Mii <c,<oo, 

where the constant Ci does not depend on z, z' G Z"^ , t G M , and e > . 

Appendix A: Proof of Lemma 6.7 

By (16.131) and fl6.15p . the function S~^^^2 can be represented as 

S;^^/,2 = ^Y1 '5l/eS + ^) ^ Ro(«£,x,l/-x) sign([ri/£] -1/1)^^/^2(^ + 1/)'^, 



where, by definition, K^^a, = e[r /e] — ex . The Parseval equality yields 
Y sign([ri/£] - yi)IU){Ke,x,y - x)glj^2{p + vY 



{2n) 



Yd 



sign([ri/£:] - ?/i)Ro(k£,^, y - x) Fy^e Gr/e^iP + y) 



dd. 



Note that 



s\gn{[ri/e\-yi)Iio{Ke^x,y-x) = {2tx) '^Fy^e sign([ri/e] -yi) *Fy^e B^o{K^,x,y-x 
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where Fy^g[Ro{Ks,x, y - x)] = e"-^Ro(K£,x, i 



9 = {9i,9) , 9 = {92, . . . ,9d) , and PV stands for the Cauchy principal part. Hence, 

• p , p ilri/e]{9i-z)+ixiz ^ 



xe''-'g'^^^,{9)*e'P-')d9 



/„i[ri/£]{ei-z)+ixiz _ _ _ . 

—— — -Ro(/s:,,,,z,^) dze'^^-'glU9re''P-')d9d9'. 



We change variables — > = — z , then denote z = 9i , and change variables 9' = 
9 — 9' . Therefore, 

S~^/,. = i^^r'^^ J e-^'^^~^>'lMPe{0) d9, (A.l) 
where Ie{9) is defined by (16.181) . and Pe{9) stands for the matrix- valued function 

/„i([ri/e]+pi)</> ^ _ 
tg(0/2) + ^ = (^1' ^")- 

Applying the partition of unity (16. 4p . (16. Sp . and formula (I6.22p . we rewrite Ps{9) in the form 

p i{[ri/e]+pi)4>p±iuJa(e-L+<t>,e)T/e^ 

^^(^) = E E y ^^^^y^^ 9(0^ + 0, ^")c^(^^i + 0, ^")* 

where cj(^)* = U^{9){I T ^C;{9))/2 . Let us fix n G M , r ^ . Since diUJ^{9) ^ for 
9 G supptyf, we can choose Eg = eg{ri,T) > such that for all e < Eg and 9 G supp(?. 

Lemma A.l Let us fix ri eM. and t ^ . Then 
(i) sup sup |Pe(^^)| < oo . 

(a) Let Ti ± diUJ(j{9)T /e ^ for all 9 G supp^f and e G (O,^^) . Then 

s 

P,{9) -2Txi^Y^ ^±iu.{e)r/e^ sign{ri ± diUJ„{9)T/e) g{9)c^{9y ^0 as e ^ +0. 
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Lemma A.l can be proved by using the technique of [T2| Lemma 8.3] or of ^ Proposition 
A. 4 (i), (ii)]. The proof of this lemma is based on the following well-known assertion: 

lim fpV [ ^^^!^^^ciz-7rze^^"(°)x(0)signa;'(0) I =0, 

A->+oo \ J Z j 

where x ^ C\ tu E , and w'(0) ^ . 

By Lemma A.l we can rewrite the r.h.s. of flA.ip in the form 

= (27r)-'^^ I e-*('-P>^/,(^)e±^^-(^)-/^' sign(ri ± d,u.ie)r/e)h^{9y 66, (A.2) 
where h^{e) = g{e)c^{e) . Substituting in f lCT) we obtain 




x(^^e±*'^-'(^)"/"%ign(ri ±9ia;,/(^)r/£)/i;^,(e)*) d9 + o{l), e ^ 0. (A.3) 

Finally, comparing the r.h.s. of ( 1A.3I) and (16.261) we see that the problem of evaluating the 
limit value of ( 1A.3P is solved by the similar way as in the step (v) of the proof of Lemma 16.61 



Appendix B: Proof of Lemma 6.8 

By (16.151) we write 

Sl/,.= J2 G'r/Ai+x)R'{^s,.,[r/e]-x,[r/e]-y)S'r/e^ip + yf, 

where ^ = e [r /e] — ea; . Change variables y ^ z = y — x and denote the series over x by 
= $e(z,r,r,/,p) 

= E ^r/e-^i^ +x)R\k,,,, [r/e] -X, [r/e] -x - z)Qli^, (p + x + zf . (B.l) 

Therefore, 

The estimate (14.91) for R and the notation (I6.14p imply the same estimate for R^ , 

|i?°(r, x,y)\< C{1 + \x- y\)-\ x,y e Z^. (B.3) 
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Next, the Cauchy-Schwartz inequality yields 

Hence, condition E6 and the estimate (1B.3P imply that |$e(-2)| < C(l + \z\)~^ . Since ■j > d , 

J2mz)\<C<oo, 

and the series in (]B.2|) converges uniformly in e (and also in T,r,l,p). Therefore, it suffices 
to prove that 

lim ^Jz) = for each z e Z'^. (B.5) 

Let us consider the series in (IB.ip . At first, note that by the property (a) for R and the 
definition (16 .Mp . the function R^{r, x, y) has the form i?°(r, x, y) = R°(r, xi, yi, x — y) , and 
R°(r, xi, 2/1, z) = for yi < . Hence, 

R°(r, [ri/e] — xi, [ri/e] — xi — zi, z) = for xi > [ri/e] — zi. 

Secondly, it follows from condition (14. 8p that W6 > 3Ks > such that for any yi > Kg 
and Vr G M'^ , |R°(r, yi, yi- Zi,z)\ < 6 . Hence, V5 > 3Ms = max{Ks, Zi) > such that 



Gr/e^(^ + a;)R°(fi;£,x, [ri/e]-xi, [ri/e]-xi - zi, z)Q^^^2{p + x + z] 



T 



xeZ'*:xi<[ri/e]-Mi 



< <5 J]|^,V(/+x)^,V(p + x + 



< C6, 



by the bound flB.4p . Let us fix 6 > . Therefore, it suffices to prove that for each l,p, z & Z'^ , 
YI 0r/e^i^ + ^)'^'i^e,x,[ri/e]-x,,[ri/e]-xi-z,,z)g'^^^,{p + x + zf -^0 (B.6) 

as e — , where := {xi G : xi E [[ri/s] — Ms, [ri/e] — Zi]} . 

At first, note that if xi E , [ri/e] —xi G [zi, M^] . Denote R^^x = ^[^/^] ■, ^ ^ R'^^"'" , 
X G Z"'^^ . Then, by property (12 .lip for Rq and property (I4.10p for R , 

sup \B^{e[ri/e\- exi,Ke,x,yi,y'i,z) -Ii^{Q,Re,x,yi,y[,z)\ ^ £ -> 0, 

for every xi G , yi.y'i G 1} , and z G U^^^ . Hence, by the bound flB.4p . we can replace 

Ke,x into (0, Ke^s) in the series (IB.6p . 

Further, we change the variable xi — > x'^ : x'^ = [ri/e] — xi in (IB.6P and denote x[ = xi , 
= (['"i/^] ~ Xi,x) . Therefore, to derive flB.6p it suffices to check that for any fixed 

l,p, z E Z'^ , and Xi G [zi, M5] , 

D, = Y G^^/^2il + Xe)B!^{0,K,^x,Xi,Xi- zi,z)g^^^^2{p + Xs + z)'^ -^0 as E 0. 
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Let d=l. Then = Ql/A^i + [ri/e] - a;i)RO(0, xi, xi - + [ri/e] - xi + z^f 

vanishes as e — ?■ . Indeed, applying the inverse Fourier transform and the decomposition 
([6221), we have 

^ — 1 J- 'J 



The last integral vanishes as e — )■ , since w^(6'i) 7^ for all Qx ^ supp (7 . 

Let d> \ . In this case, we put, for simphcity, R°(0, Kg^^j, • . •) = R^lO, /^e.s, 2:1, Xi — ^i, . 
Using the Fourier transform we rewrite as 

Let us change variables 9 = {9i,9) , 9' = {9[,9') as follows: 9' —> 1^ = 9 — 9' , and denote 
6'i = v?]^ , 6*^ = 6*1 . Hence, 

D, = {27i)-'^+^ J e-'^^-P~~'>^A,{9)B,{9)d9, (B.7) 

where 



Ai0) = 


2vr ^ 






(B.8) 


Bs{e) = 


2vr ^ 




eyd^u 


(B.9) 


UO) = 


{2n] 


Jd-1 


J2 e'^m°(0,£[f/£] -ex,.. 


.))(i<^ (B.IO) 



(cf. formula (I6.18P ). Below we prove the following bound for le{0) (cf- 06.24P ): 

s 

^^(^) = E E (^) e^'-(^)^/^' A^<,(r, f; 9) + o(l), e ^ 0, (B.ll) 



<T=1 ± 

where h^{9) is defined in (16.221) . 



Al,{T,f;9)= j RO(O,fTV,-a;,(0)r/£-y,...)/C^(r,y,0)rf^, (B.12) 
/C^(r, I/, 0) , r > , 1/ G M'^"^ , 9 eT'^ \ C , stands for the following matrix-valued function 



g±j7rs/4 gTi/(2T)y-{V2a;<,(e))-is 
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and s denotes the signature of the matrix (V|a;^(^)) , 6 = {01,6) e T'^ \ C (cf. ^Ml). 
To prove ( IB.lip we first apply the Poisson summation formula (see f l6.19p ) and obtain 

J2 e^^-^R°(0, e[f/e] - ex, . . .) = e^'^^^ Yl "'^/^ ~ 27^/6, . . .)e'^'/'^-^, 

where R°(0, . . .) stands for the Fourier transform R°(0, . . .) = Fj;^(p[R°(0, f, ...)]) 
f,^e M*^-^ . Conditions fim]) and fim]) yield 

R°(0, - 27rfi/e, . . .) ^ as e ^ (B.13) 

uniformly in <^ G [— vr, vr]'^"^ . Hence, by f lB.10|) and f IB.lSp . 

Ue) = {2718)-'+' j e^^-('+[^>l)^f/^,(^i,^-^)R°(0,-^/£,...)# + o(l) 

= (27r)-'^+i j e-*^^-('+[^>l)^f/^,(^^i,^ + £(^)R°(0,<^,...)ci<^ + o(l). 

[-7r/e,7r/e]<^-l 

Secondly, using the representation (16.221) for Qf{0) we rewrite 1-^{0) in the form 

X,{e) = (27r)-'^+^^ j e-'^^-('+[*'/^])e^'"'^('^''"+"^)"/^'/i;^(^i,^ + £(^)R°(O,(^,..0ci^ 

+o(l) = (27r)-'^+i ^ /i;^(0) y e-^^-^V^"'^(^i'^'+=^^^/^'R°(0, . . .) # + o(l), 

by conditions (I2.1ip and (I4.1ip . Finally, we use the Taylor sum representation for Ua{9i,9 + 
Eip) (cf (16.231) ) and obtain the bound ( IB. Ill) , since 



if,...) dip 



= ^,--i.^v,..(.)./JR°(0, ^, . . .)e±^(-/^)^-(^i-(^))^] = AUr, r; 9), 

where Af^{T,f; 9) is defined in (]B.12p . The bound (IB. lip is proved. 

Now we prove that the r.h.s. of (IB.7P vanishes as e , using the Lebesgue dominated 
convergence theorem. At first, we show that 

1^,(^)5,(^)1 < v{9), ye > 0, where v{9) G L\T'^-'^). 

Indeed, by conditions 14' and (14. lip . 

sup sup |^^^(r, f;9)\<C<oo. 

e>o eeT^ 
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Since h^{e) = g{e)U^{e){I T iCa{0)) /2 , the bound flRTTjl yields 

s 

\Iem <Y.^Ciu;\9) + C2U^{9)) + Cs, 9er\C. 

(T=l 

Write v{e) := / (I + \\V-\e)\\) dOi , 6 E T'^-^ . Therefore, by dEl and (iROl) . 



1^(^)1 < C j \Iemd9i<Ciiv{9)Y/^, 

XI 

\Bem < c j \gi^^,[e)\der<cMW^ 

XI 

and, evidently, \Ae{e)Be{9)\ < Cv{e), where v{e) G ^^(T"'-^) by condition E6. Therefore, 
to prove that the integral in flB.7p tends to zero it suffices to check that -Be(^) vanishes as 
e — 7- for a.a. 9 E f^^i . By the representation fl6.22p . 

2^ „i I J 



For fixed ri G M"*^ , r 7^ , and G 1} , mes G : (9i(X'o-(v9i, ^) = 0} = for 

a.a. fixed 9 E T"'^^ . Hence, B^{9) vanishes as e — )■ . It can be proved similarly to the 
Lebesgue-Riemann theorem. ■ 



Appendix C: Proof of Lemma 6.6 in the general case 

We prove Lemma [6.61 under the weaker condition fl2.13p - fl2.14p on Ro(^)^) than (12. lip . 

Step (i) of the proof is not changed and we derive formulas fl6.16p - fl6.18p . 
Next step is to apply the Poisson summation formula f ]6.19p . However, in general, we can 
not apply fl6.19p to Ilo{r,9) . Therefore, introduce 

RrAy^(^) = ^oiy,0)e-''^y-'^\ Z/GM^ 9eT', (C.l) 

with P E {d + 2,2{6—l — d)/d) and 6 from fl2.14p . The asymptotics of fl6.18l) is not changed 
if we replace the function Ro(?/,6') by R.,.£(?/,6') in fl6.18l) . Indeed, consider the following 
series: 

J2 J e*('+^)-^^f/^,(e - ^) d^(Ro{e[r/e] - ex, 9) - R.,,(£[r/e] - ex, 9)^) 
= E •••+ E ■■■ = M0) + Bs{9). 
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and prove that the series ^4^(6') and Bs{0) vanish as £ — ?■ . For the first series As{6) , we 
apply Lemma [6.41 (i), definition (IC.ip and property II and obtain 



\x\<CT/e^ 

uniformly in 6 , since f3 > 2 + d . For the second series Be{6) , Lemma [6.41 (ii) yields 



\x\>CT/e'^ 



X] 



0, e ^Q. 



Step (ii): Now we apply formula (16J9D to ^ e"-^Rr,e{4r /e]-ex, 6) . By RrA^, 9), s G R"^ , 
6* G T*^ , we denote the Fourier transform of Rr,e(y,6') w.r.t. y (see (12.121) ) and obtain 

= (27r£)-^ f e*('+['-/-])-'.^^5^^^ ^0_^)(^J2 ^rA-v/^ - 2™/^, 0)) + o(l) (C.2) 



(cf (I6.20p ). We check that the contribution of the series with n 7^ in flC.2p vanishes, 



Rr,ei-^/^ - 2vrn/£, ^) ^ as e ^ 0, 



(C.3) 



uniformly in G T"' , ip E [— vr, ttY (cf (I6.2ip ). At first, we derive the formula for Rr,£(.^, 



Rr,.(e, 0) = e^"-«(47rA;)-^/2 j ^-is-r^-\s-e/m^(^Q^ ds), 6 G 



(C.4) 



where k := and ^{9,ds) from (I2.13p . Indeed, by definition (IC.ip . 

R^_,(^,0) = y" e*«-^-'^(^-'-)'Ro(i/,^) rfy = e"-«-l«l'/(^^) J e-''^y~'-'^/^^^^^'Ro{y,0)dy. 

On the other hand, by definition (I2.13p . for any a G M'^ , 

J e-'=(^-")'Ro(?/, 9) dy = {2Tr)-'^ J J e-'y'e"^^^-"^' fx{e, ds) dy 



(2vr) 
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Hence, 



oir-i-\i\vm 



d/2 / is.(r+i5/(2fc))^-|s|2/(4fc) 



/x(6', ds). 



This implies formula (IC.4I) . 

Now we prove fIC.Sp . Let us set ^n,e = — 27™/^, where < vr , i = 

apply flC.4p and devide the integration into two: \s\ < TT/{2e) and \s\ > it/{26) : 



= /1 + /2, 



(C.5) 



where 



|s|<7r/(2e) 



(C.6) 



|.|>7r/(2£) 



Let |s| < 7i/{2e) . Put = -^Pi/e - 2'nni/e . Then \si-in„e\ > \27i\ni\/e - \si + ipi/e\ \ > 
{7c/6){2\ni\ —3/2) if rij 7^ . Therefore, there exist positive constants a,b such that 



Hence, there exist positive constants and bs such that 



s=0 



since k = e^^ with /3 > . Let |s| > 7r/(2£:) . Then 

e-{s-^n,e)Vm ^ g-(^V(fc£2))(n+^/(27r)+.£/(2^))2 ^ C'^(ev^)'^ + (:72- 



(C.7) 



Using dUj]), ( lUTl) . and condition (EH]), we obtain 

h ~ (Ci + C2ieVk)-'^) sup / < C'3e^-i-(^/2+i)d ^ ^ _^ 

|s|>7r/(2£) 

because k = and 5 > 1 + (/3/2 + This completes the proof of the bound (IC.Sp . 
Step (in): We apply flC.2p . fIC.Sp . change variables (p — j- —eip and use formula (16.221) : 



7,(0) = (27r£)-^ J e*^-('+[^'/^l)^f/^,(e-<^)R,,,(-<^/£,^)d<^ + o(l) 
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(27r)-^ j e-*^^-('+[^/^])^f/^,(^ + ev)RrAv, 9) + o(l) 

[-7r/£,7r/e]'* 

{2tt)-'^ J2 J e-*^^-('+[''/^l)e±'"'^(^+^^)"/^'/i^(^ + eip)RrAv, 6) dip + o(l). (C.8) 



The asymptotics of (IC.SP is not changed if we replace e *^'^ ('+[''/^]) by e in the integrand. 
Indeed, we check that the integral 

D,{d) = j (e-^^^-('+[^/^l) - e-i^-r^e±iu^Ae+e^)r/e^^T^g + eip)RrA^, 9) dip 

[—n/e,TT/e]'^ 



vanishes as e — )■ . We apply formula (IC.4p and change the order of the integration 
fi{9, ds) dip — )■ dip fi{9, ds) : 



[-vr/e.Tr/e]'' 



""((4^^ y e-^^-^e-l^-^l'/(^^)/x(^,rfs)jrf¥. = j e-^^-^C,is,9)fi{9,ds),{C.9) 

where Ce(s,^^) stands for the inner integral in (iy) : 

C,(s, 9) = J e-l^-^l'/(^'=)(e'^-('-"[^'/^l-^') - l)e^''^-^'+'^^^/'' h^{9 + eip) dip. 

[~iT/e,iT/e]'^ 



Note that 



\Ce{s,9)\<esnpm9)\—^^ [ e-\^--\'/^''^\ip\dip, 



where 

1 



[— Tr/SjTr/e]'* 



Hence, by (JCS]) and fimj) . 

supegTd l-Oe(^)l — C'e. Similarly, the asymptotics of fIC.SP is not 
changed if we replace h'^ij) + Sip) by h^{9) in the integrand. Therefore, 

h{9) = Y,K{d)NZAr,r;9) + o{l), 5^0, (C.IO) 

where N^^{T,r; 9) stands for the matrix- valued function, 

N^^,{T,r;9) = {2tt)-'^ J e-''^-''e^'''-^^+'^^^'''Rr,e{v.O) d^. (C.ll) 

[— 7r/£,7r/e]'^ 
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Step (iv): Now we apply the Taylor sum representation (cf (16.231) ) and replace ujcj{0 + eip) 
by uj„{0) + eVuj„{0) ■ ^ + ■ H„{0)^ . To do it we check that 



-7r/£,7r/e]° 



as e — )■ . As in Step (in), we use formula (IC.4P and change the order of the integration 
fi{9, ds) dip —J- d(p fi{9, ds) . Therefore, 

L,{e)= I e-''-'-D,is,9)ii{9,ds), 



where Di.[s,9) stands for the inner integral in dip 
1 



DJs, 



dip. 



Moreover, by condition (12.141) . 



sup 



e-'''D,{s,9)fi{9,ds) 



\>N 



\s\>N 



'^\'/^^''^dip^\^i\{9,ds) 



< C{1 + N)-^+^ ^0 asAT^oo, 
uniformly in e . Hence, it suffices to prove that for any fixed G N , 

L^{9):= j D,{s, 9) fi{9,ds) ^0 as e -> 0. 

\s\<N 

To check ( ]C.12|) we estimate the inner integral -0^(5, 9) : 

1 f _|„_,„|2, 



(0.12) 



ip\ dip 



-n/e,iT/e]'^ 



Hence, for fixed A^ , 



\L^{9)\<Ce 



:^''/^+\sfM{9,ds)<C2iN)e, 



\s\<N 



and (1C.12|) follows. Therefore, sup^g-ipd |L£(6')| — )■ as e — )■ , and 



(C.13) 
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where, by definition, 

= jRrArTVuj^{0)T/e-x,e)K^{T,x,e)dx, (C.14) 

with the function K^{T,x,e) ( r > , x eR'^, G T"' \ C ) from flO^ . 
Remark. The integral N^^{T,r]6) can be rewritten in the another form. Namely, 

N^^Ar,r;e) = e^'--^'^^/^'C^^^ir,r;e) + oil), e ^ 0, (C.15) 
where, by definition, 

Cf^,{T,r;e) := {2n)-'^ [ e-^''<'^^^-^^^^/'^ e^'^^^''^^^^'^^' ^{6, ds). (C.16) 



Proof. We first substitute formula flC.4p in fIC.lip and change the order of the integration 
fi{9, ds) dip —J- dip fi{9, ds) . Therefore, 

N^^Ar,r;e) = (2n)-'' I e-- [j^^ J e-^'-^^'/^^'^e^^-^^'^^^^^/^' dip) f^iO, ds). (C.17) 

Rt* [-7r/e,7r/e]<^ 

Secondly, we replace ijJa-{9 + eip) by ijj^{9 + es) in the inner integral. To do this we devide the 
integration into two: |s — > e'^ and |s — y^l < e'^ with a 7 , 7 G (1, (3/2) . For |s — y?! > , 

e-'^-l^/^^'=^rf^<^7-r^ - <C^^-.:t7^^0 (C.18) 



as 



[-7r/e,7r/e]<^,|s-v?|>e^ 

£ — )■ +0 , since k = and 7 < /3/2 . Note that 



1 



(47rA;)^/2 

[— 7r/e,7r/e]'' 



e-'^-^''/(^'=)c/(/^ = 1 + 0(1), £^0, 



where k = . Therefore, for \s — ip\ < , we estimate the difference 
1 



{ATikyi^ 



[-7r/e,7r/e]'*,|s--(/3|<eT 

^ 7T-tW / e-l^-^l'/^^'^^^^+V/e^ c/^ < Ce^-\l + o(l)), (C.19) 



(47rA;)^/2 

[-7r/e,7r/e]'* 
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where the last expression vanishes as £ — )■ since 7 > 1 . By flC.lSP and (IC.19P , the inner 
integral in flC.17p is 

1 



[-7r/e,7r/e]'*,|s-(/3|<eT 

= / + 0(1) 

[-7r/e,7r/e]'*,|s-v5|<eT 
[~n/e,Tr/e]'^ 

uniformly in 6 ^T'^ . We substitute the last expression in fIC.lTp and get 



N^^^{T,r;e) = (271)-'^ J e~''-'e^'''^^'^-^"^^/'' ^{6, ds) + 0(1), e ^ 0, 

uniformly in 6 E T'^ . Finally, we apply the Taylor sum representation f l6.23p to uJcr{6 + £s) 
and obtain ( ]C.15I) -( IC.16I) . The integral in ( ]C.16p can be taken over |s| < with some 



G N, by condition iHA^i . 
Step (v): We substitute /^(^) of the form (IClTnl) with Nf^„{T,r;e) from (1(113!) (or (EH])) 
in the r.h.s. of (I6.17p . Applying the decomposition (I6.22p to ^^y^2(^)* we obtain (I6.26p . with 
BJ^^ (or Cf^^ , respectively) instead of . It remains to study the behaviour (as e — )■ ) 
of integrals of the form fl6.27p with (or , resp.) instead of . 

Let oJaiO) ± ^a'{(^) ^ const± . In this case, the oscillatory integrals 

I^^,{e) = (27r)-^i J e-'('-P)-^e*("''(^)±"-'(^))"/^'/i;(0)C-,(r,r; e)h^,{e)* dO (C.20) 

vanish as , since sup |C~^(r,r; 6')| <C < 00 , and hence, h~{6)C~^{T,r;6)h^,{6)* E 

L^iT'^) . The identities u:rj{9) ± ujfji{9) = const-t in the exponent of (]C.20p with const-t 7^ 
are impossible by condition E5. The identity uja{0) +uja'{0) = implies uja{0) = uja'{0) = 
which is impossible by E4. Therefore, if a ^ a' , /^^, (e) = o(l) as e — )■ . If a = ci' , 
I^„{e) = 0(1) , and only I^^{e) contrubutes in the limit. Hence, (cf (16.280 ) 

^;./.. = W-'l J e-^^^~^y'J2h^{e)BXAr,r;e)h^{erde + o{l), e ^ 0. (C.21) 



J. 



•d 



Step (vi): To obtain (I6.28p . it remains to replace Bf^ by Af^ in (]C.2ip (or R^^ by Rq 
in (]C.14p ). We check that the difference B~^ — vanishes as e — i- (for the difference 



^ta - Ka the proof is similar). Indeed, by (l37[3D . dOTl) . (IClij) . and ( lOSj) 
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with C{9) = e^"^/^(27rr)-"'/2|det/J^(^)|-i/2_ ^^ite c = 2|r| max |Va;^(^)|. We devide 

the integration into two: \x + 'VijJa{0)T/e\ < c/e and \x + 'Vuja{0)T/e\ > c/e. For \x + 
Vuja{d)T/e\ < c/e, 



\x+^Ua(B)T/e\<c/e 



< c 



e"\x + Vuj^{e)T/e\' dx < C^e^-'^^'^ ^0, e ^ 



|a;+Va;<^(e)T/£|<c/£ 

since /3 > + 2 . For |x + Vco'o-(^^)r/£:| > cje , we apply the integration by parts in every 
variable Xj , i = l,...,d, For simplicity, let d = 1 and /i stand for —1/{2t)H~^{9). 
Therefore, 



+ 00 



c/£-<(e)r/£ 

1 

2ih 



Ro(r - u;;(e)r/£ - x, ^) (l - e-^''(-^W-/^+-)')e'^"' 



-Ro(r-c/£,^) 1 - e 



a-={c-<{e)r)/£ 



a,Ro(r - cj;(^)r/£ - x, 9) (^1 - e-^''(-^W-/^+-)') ^ 



dx 



Ro(r - uj'^{e)T/e-x, e){uj'^{e)T / e + x)e-^'("^(^)"/^+")'^ rfx 



X 



(c-<(e)r)/£ 
+00 



+ / Ro(r-a;;(e)r/£-x,^)e-"'("^(')"/^+")'^dx 



X^ 



(c~<^;(e)r)/e 

= h + h + h + h- 

The terms Ji and /4 tend to zero as e — )■ , since |Ro| < C < 00 by condition II. Note 
that |x| >c/{2e) if \x + uj'^{9)t / e\ >c/e. Hence, 



2 1 



/3I < Ci£^ / {uj'Mr/e + x)e-''^^''^^'^^/'+^'>"^^dx 



X 



(c-a;;{e)r)/e 
+00 



< C^e^^^ j ye-''y"dy = C, e e-^''(^/^)' ^0, e ^ 0. 

c/e 

To prove that the contribution of the integral I2 vanishes, we repeat the integration by parts 
and use the bound sup |9^Ro(r, ^)| < C < cxd , where k G [0, d'^/2 + 2d] . This bound 



follows from condition fl213il - fl2lD . 
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Appendix D: Local conservation law 

Let v{x, t) be a solution of (12. ip with the finite energy. The local energy in the point x G Z*^ 
is defined as 

We derive formally the expression for the energy current of the finite energy solutions f (x, t) . 
For the half-space fi^ := {a; G Z"' : x = (xi, . . . , x^), a;^ < 0} , /c = 1, . . . , d , we define the 
energy in the region VL^ as 

By formal calculation, using Eqn (12. ip we obtain En^if) = — ^ jk{x',t), where 

jk{x',t) stands for the energy current density in the direction = {6ik, ■ ■ ■ ,6dk) , k = 
x' eZ"^ with x'fc = 0, 



jk{x',t) := ^ { X] vix' + mek,t) -Vix' + mek-y' -pek)v{y' + pek,t) 

y'eZ'^-y'k=o m>i,p<o 

- ^ + mefc,t) ■ ^(x' + mcfc - -pefc)?;(?/' + pefc,t)|. (D.l) 

m<0,p>l 

Now let v{x,t) be the random solution to (12. ip with the initial measure /Iq satisfying 
VI and V2. Therefore, for any x e Z-^ , r G R \ , and r e M'^ , the average energy is 

Eg[^(x + [r/e],r/5)] = K^ix + [r/E],0)] 

= ^ tr [Q^,V/,([r/5] + X, [r/e] + x) + QZ/ei[r/^] + x, [r/e] + y) ■ V^{x - y) , 

by condition El and the uniform bound for the correlation functions of the measure /x^^^ 
(see [HI Lemma 5.1] or Lemma [7. ID : 

sup sup sup IIQJ^ , (;z,z')ll < < oo- (D-2) 

e>0 jj=0,l 

It follows from Theorem 13.21 that for any r 7^ , r G M'^ , Eq[£^(x + [r/e], r/e)] — t- e(r, r) as 
e — > , where 

e(r,r) = ^ tr [g^;,(0) + ^ (x)\/^(x) 

= \{2^)-'ir J + de = (27r)-^tr J ^(0) dO. (D.3) 

The last equality follows from condition E2 and Remarks 13.41 (i). 
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Similarly, by flD.l|) we obtain 

y'gZ'*:yJ,=0 rn<0,p>l 

- ^ tr + mefc,|/' + pefe) ■ - y' + (m-p)efc 



Therefore, by Theorem 13.21 the following limit exists, limEQh'A;(x+ [r /e],T/e)] = jA:(T, t) for 
any t , r eW'- , k = 1, . . . , d . Here 



jk{.'r,r) 



1 



^ ^ ti^ql^,{x'-y' + {m-p)ek)V^{x'-y'+{m-p)ek) 

y'eZ'i-y'k=o m<-i,p>o 

^ tr (x' -y' + {m- p)ek)V^{x' - y' + {m - p)ek) ) • 



m>0, p< — 1 



Write x' — y' =: z' G Z'^ with z'f. = , m — p := s E ^ and change the order of the 
summation in the series. Therefore, 



(D.4) 



--(27r)-^tr / qll{o)dkV{e)de, k = l,...,d. 



Write j(r,r) = (ji(r,r), . . . ,jd(r,r)) . Finally, flD^l) . flOil) and Corollary |3J] yield 

a^e(r,r) + V,. ■ j(r,r) = 0, r G M, r G M'^. 

Remark. In Section 2.3, we give the example of the "local equilibrium" initial measures 
fif) . Namely, let qQ (z) be the correlation functions of the Gibbs measure g (see Definition 
ESD with /3 = 1, i.e., = V-\e), ql\e) = J, qf{e) = ql^{e) = 0. Put Ro(r,^) = 

T{r)qo{9) , where the function T{r) , r G M'^ , is defined in Section 2.3. Moreover, let /Iq , 
£ > , be Gaussian measures with the correlation functions Ql^{z, z') defined in fl2.17p . Then 
conditions VI and V2 hold. In this case, the limit correlation matrices $^^{9) have a form 

(T=l 

S 

cr=l 

where, by definition, T^(r, r; 9) := (l/2)(T(r + Va;^(^)r) ± T(r - Va;^(e)r)) . Therefore, 
e(r,r) = (27r)-'^^ /" T^(r, r; tr n.(^^) d^. 
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In |6], the "locally conserved' quantities were studied in the case when d = 1 . Let us 
consider these quantities in many-dimensional case. At first, introduce the following matrix- 
valued functions 



E{z + h,z;Xo) := ^(^Vi{z + h) ^ Vi{z) + Vo{z + h) ® ^ V{z - z')vo{z'] 
A{z + h,z-Xo) := - {vQ{z + h) ®Vi{z) - Vi{z + h) ® Vo{z)^ , z,he Z'^, Xq = {vo,Vi). 

The "locally conserved" quantities Xl{ip, Xq) , Y^{ip,Xo) are defined as follows. For ip G 

C^{R'^), we set 

Xli^,Xo) := e''J2^{ez)Eiz + h,z;Xo), 



Theorem D.l Let conditions 11-14, VI, V2 hold. Then for any G Cl{W^) , /i G , 
r 7^ , there exist the limits 

\imK/e[Xl{v. ■)] = E(^; r, h), limE^/jy,^(¥,, ■)] = A(^; r, h). 

r/ie matrices E(v9; r, /i) and A((y9; r, /i) are g'Zi'en 6y their Fourier transform in the following 
way. Write 

E(y.; r, e) := e*^-^E(y.; r, h), A(^; r, ^) := e*^-^A(y.; r, /i), ^ G T'^. 

Then 

E(y,; r,e) = ^{r) (g^^^) + = / vir)^^ dr, (D.5) 

Mcp; r,e) = y v^(r) (g^^^) - ql^M) dr = j ^ir)q'^M dr, (D.6) 

where ql^r{9) are defined in liS. 

The proof of Theorem D.l is based on Theorem 13.21 and the bound (lD.2p . 
Remark It follows from f lD.Sp . flD.6p and Corollary 13.31 that E and A satisfy the equations 

dM^; r, 6) = lio^ie) J Vio^ie) ■ Vip{r)f^M dr = too ^{e)V 00^(6) ■ A{V^; r, 6), 
drA{ip; T, 9) = ~iuj-\9) j Vuj,{e) ■ Vv{r)ql]M dr = -iuj-\e)VujM ' E(V</^; r, 9). 
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